Quasi exactly solvable (QES) equations with 
multiple algebraisations * 



Yves Brihaye 

Faculte des Sciences, Universite de Mons-Hainaut, 
B-7000 Mons, Belgium 

Betti Hartmann 

Department of Mathematical Sciences , University of Durham 
Durham DH1 3LE , United Kingdom 

February 1, 2008 



Abstract 

We review three examples of quasi exactly solvable (QES) Hamitonians 
which possess multiple algebraisations. This includes the most prominent ex- 
ample, the Lame equation, as well as recently studied many-body Hamiltoni- 
ans with Weierstrass interaction potential and finally, a 2 x 2 coupled channel 
Hamiltonian. 



1 Introduction 

Quasi exactly solvable (QES) operators can be formed by considering elements of the 
enveloping algebra of some Lie Algebras realised in terms of differential operators 
preserving a finite dimensional vector space of smooth functions. Some of them can 
be transformed into Schrodinger operators after performing a suitable change of vari- 
able and a suitable change of function (also called a "gauge transformation"). These 
Schrodinger operators admit a finite number of algebraic eigenvectors and are thus 
called QES. For certain QES Schrodinger operators there exist several choices of the 
"gauge functions" which then lead to different algebraisation of the corresponding 
eigenvalue equation. If this is possible, we speak of "multiple algebraisations". The 
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generic example of such an equation is given by the Lame equation. In this contri- 
bution, we present further examples of QES equations which have multiple algebrai- 
sations. These are equations describing N bodies interacting through a Weierstrass 
potential and a system of coupled QES equations related to the sphaleron stability 
in the Abelian-Higgs model. 

2 Quasi exact solvability (QES) 

On the one hand, we have a Hamiltonian H(x) = — + V(x) with x e M. d and the 
corresponding eigenvalue equation Hip{x) = Eip(x). On the other hand, we have a 
Lie Algebra Q of dimension rig realised in terms of operators Jk = JkijJ, V^), y e ~R d , 
k = 1, ...,ng which preserves a m-dimensional vector space V of smooth functions. 

Now, if there exists a) a change of variables x = x(y), b) an element of 
the enveloping algebra of Q as well as c) an invertible scalar function fix) such that 
H{x) = /(a ? )7i(y)/~ 1 (a ? )|^ = ^), then the Hamiltonian H preserves the vector space 
fV: HfV C fV. As a consequence, the restriction of the operator H to the vector 
space V leads to an algebraic equation [1]. This equation allows to construct m 
eigenvectors algebraically. In this construction, the factor / which is called "gauge 
factor" plays a crucial role. The above procedure is usually called an algebraisation 
of the Schrodinger problem. 

In this contribution, we will put the main emphasis on quantum Hamiltonians 
H which possess multiple, i.e. several algebraisations. More specifically, more than 
one "gauge factor" /, say f a , a = 1,2,3... exists. This allows to write H(x) = 
f a (x)%(y)fa 1 {x)\y=y(x) where Jk >a are now the generators of the Lie Algebra Q a . Cor- 
respondingly, the invariant vector spaces V a of the generators J^^ can also be of 
different nature. 



3 Hamiltonians with multiple algebraisations 
3.1 The Lame equation 

A natural example of a QES equation for which multiple algebraisations exists is the 
Lame equation: 

(~C^2 + N(yN + 1 ) A;2sn2 ( ;r ' = E ^ X ) W 

where sn(x, k) is the Jacobi elliptic function of period AK(k) with K{k) being the 
complete elliptic integral of first type. The potential in this problem is periodic with 
period 2K(k) and the parameter N determines the height of the potential. For the 
equation (1), four different algebraisations (a = 1,2,3,4) exist. For the case of even 
iV = 2n, neN, the gauge factors are given by: 

/i — 1 , fi — sn(x, k)cn(x, k) , / 3 = sn(x, k)dn(x, k) , / 4 = dn(x, k)cn(x, k) . (2) 
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The relevant change of variable is y(x) = sn 2 (x,k). Then H a = f a Hf a 1 can be 
written in terms of the operators 

2 d d n a d 

J + = y Ty- nay > J « = y Ty-^' J - = Ty> (3) 

where n\ = n, n 2 = = = n — 1 and Vi = V(n), V 2 = V3 = V4 = V(n — 1). 
Here V(n) denotes the set of polynomials of degree less or equal to n in y. The four 
underlying algebras Q a are in this case equal to si (2). 

3.2 Many-body Hamiltonians 

In this section, we will discuss how the idea of multiple algebraisation can be applied 
to many-body Hamiltoninas with N degrees of freedom. The Hamiltonians considered 
here are generalisations of Olshanetsky-Perelomov- [2] and Inozemtsev [3] -type. More 
specifically, the Hamiltonian is given by: 

N 

H N (x) = -A N + a(a - 1) ^ (P(xj + x k ) + P(xj - x k )) 

i,k=i 

N N 

+ Ab(b-l)J2 P ( 2x k) + cJ2 P ( x k + iP) ■ ( 4 ) 

k=l k=l 

P( z , 92-, 93) denotes the Weierstrass function of parameters g% which are related 
to the periods on the real and imaginary axes of the complex plane, 2a and 2i/3, 
respectively, a, b, c are coupling constants. For later use, we recall the following 

identities: (i) (^f^j = 4P 3 (z) - g 2 P(z) - g 3 and (ii) P(x + i(3) = e 3 + (e 2 - 
e^)sn 2 {^/ei — e^x, k) with k = f^Eff- The singularity of the Weierstrass function 
at the origin ( P(e « 1) = ^ + 0(e 2 )) renders the potential singular for \x k — 
x\\ — > such that the eigenvalue problem can be restricted to the domain D = 
{(xi,x 2 , ■■■Xn)\ < X\ < x 2 < ... < xn < a}. The condition that the wave function 
vanishes on the boundary of D leads to a discrete spectrum of H^. 

Before addressing the main result of [4], we give the following definition of sym- 
metric polynomials [5]: 

T\ = Zi + Z 2 + ... + Z N , T 2 = Z\Z 2 + ZiZ 3 + ... + ZN-lZN , •••• , T/V = ZiZ 2 Z 3 ...Z N (5) 

and the definition of the vector space 

N 

V m = span{T[W l 2 \...T l » ,Y, l i<™}, dim{V m ) = C^ +m (6) 

where C denotes the combinatoric symbol. It is well known [5] that the space V m is 
preserved by a set of operators realising the Lie Algebra sl(N + 1). 
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Table 1: Specific limits of the five parameter Hamiltonian H(a, b, c, e 2 , 63) 



prameters 


# of algebraic eigenvectors 


Remark 


(a,b,c, e 2 ,e 3 ) 







(a,fr,c m ,e 2 ,e 3 ) 


AT 


see [4] 


(a,0,c m ,e 2 ,e 3 ) 


/~iN+m 1 r,y~<N— 1+m 
Cyjy -t-OCyjy 


see further text 


(0,0,c m ,e 2 ,e 3 ) 


(4m + l) N 


N decoupled Lame equations 


(0,0,c m ,e 2 ,e 2 ) 


infinite 


N free oscillators 



Proposition 1 : Let 

N N b 

(i) fi(x) = H [P( Xj + i(3) - P{x k + i(3)} a H [P'(x k + i/3)] 

j<k k=l 

(ii) z k = P(x k + i(3) 

(iii) H N (z) = ii~ 1 (x^H n (x)ii(x)\ £=£{s) . 

(i) is the "gauge factor" and (ii) the new variable. Then Hn(z) depends only on the 
symmetric variables t±,..., t n and H-N^Vm C V m iff 

c = c m =[2m + 2a(N - 1) + 46] • [2m + 1 + 2a(A^ - 1) + 2b] , me N . (7) 

If in the above expression for the coupling constant c the parameter m is an integer 
the total number of algebraically obtainable eigenvectors is C^ +rn . 

The five parameter Hamiltonian H(a, b, c, e 2 , e 3 ) has well-known limits for specific 
choices of the parameters as is denoted in Table 1. 

It should be remarked that out of the (4m + l) Ar eigenvectors for the case a = b = 
a total number of C^ +m + 3C^~ 1+m are completely symmetric under permutations of 
the coordinates and as such can be written in terms of the r-variables. This number 
contrasts with the number of algebraic states existing in [5]. This gives a hint to 
the fact that further algebraisations might be possible to construct by introducing 
different "gauge factors" . 

In [6], we introduced these additional factors of the form 

N 

fc=i 

such that Hn(z) = j J '~ 1 {z)'^N{.z)ji{z). This then leads to the following 
Proposition 2: If the exponents u i: i — 1,2,3 in (8) are chosen according to Vi = 
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or z/j = | — b, then li.N{z)M-m Q -Mm where m :— m + (b — \)rif. 

In the above definition of m, rif denotes the number of non-zero values chosen for 
the exponents i^, that is to say rif = or 1 or 2 or 3. This proposition results into 
the existence of seven possible new algebraisations. The algebraisations studied in 
[5] correspond to the choice rif = 0. The choices rif = 1 or rif = 2 lead to three 
new algebraisations according to the possible permutation of the non-trivial factor. 
Finally, the choice rif = 3 obviously leads to only one new algebraisation. In order to 
avoid singularities of all possible "gauge factors", we will restrict b to < b < \. 

It should be remarked that the previous condition me N is now replace by the 
condition me N. Then, for fixed b and rif the parameter m defining the coupling 
constant c has to be chosen appropriately. Though m is still quantised, the number 
of possible values of the coupling constant c = c m is considerably increased. To 
finish this section, we investigate for which fixed values of the coupling constants 
(a,b,c,e2,es) multiple algebraisations exist. For generic values of b one (resp. three) 
algebraisations occur for the case when m is an integer and ny = 0orn/ = 3 (resp. 
rif = 2 or 3). For 6 = four algebraisations coexist corresponding either to (a) rif = 
(with invariant subspace A4 m ) and rif = 2 (with invariant subspace M. m -\) or to (b) 
rif = 1 and rif = 3. Finally the case 6=1/6 leads to two algebraisation with rif — 
and rif = 3 with respective invariant spaces M. m and M. m ^\. 



3.3 Coupled channel equation 

Apart from studying many body equations with algebraic properties also systems of 
coupled Schrodinger equations of the form 

H ~ ^ II + ( Vu{x) Vu{x) \ (9) 
H2 * 2 -~^ 2+ \V^x) V 22 (x) ) (9) 

can be studied. In [7] a classification of equations of this type which admit multiple 
algebraisations was achieved. Here, we will study one example which has a physical 
motivation. The potential reads: 

V = sn 2 (x,k) ( A ~ 2b A l 2 b) +b (o ® 1 ^+ dsn (x,k)cn(x,k)(^° 1 J) (10) 

where k e [0 : 1], b e R, me N and A = f (4m 2 + 2m + 1), 9 = 4b 2 - k 4 (4m + l) 2 . 
Using again the change of variable y = sn 2 (x, k) and an appropriate "gauge factor" 
which is now a 2 x 2 matrix, we were able to show that the Hamiltonian H 2X 2 has 
a quadruple algebraisation. Here we present two of the algebraisations and refer the 
reader for the remaining ones to [7]: 

(F-iHF) { c { V(m-l)\ _ ( sn 0\/l Kl \ 

{ J V V ( m ) J V V ( m ) J \0 cn)\0 l) 
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(G HG \v{m-l) ) ~ \V{m-l) )' & - \ dn){o 1 ) 

where K\ = —0/(2b+k 2 (l+4m)). Note that in suitable realisations (which we will not 
present here) the operator F^HF can be interpreted as an element of the Super-Lie 
Algebra osp(2,2), while G~ X HG is an element of the Lie Algebra sl(2) x sl(2). This 
provides an example of multiple algebraisations where the different algebraisations 
are related to distinct Lie Algebras. 

The limit k — of these coupled equations is non-trivial and leads to the following 
operator: 

it d2 it , oh ( cos 2 (x) - | cos(:r) sin(x) \ r [ft „, r , ni x 

i/2x2 = -^ ll2 + 26 Uos( a ;)sin(J) sin 2 (i)-i ) ' ^ ^ [0, 2JVtt] . (11) 

In this limit, the dependence of the potential on the integer m disappears and the 
system becomes exactly solvable. In order to construct the flag of invariant vector 
spaces, it is useful to introduce the vectors 

It is easy to show that the following vector spaces are left invariant by the operator 
(11): 

Vb = Span{£(iV)} , V = Span{G(iV)} 

V k = Span{£(iV + k), E(N — k)} , fc = 1,2,3... 

V k = Span{G(A^ + k), G(N — k)} , k = 1,2, 3... (13) 

The operator can be diagonalized on any of these vector spaces and leads to the 
following sets of eigenvalues. Using natural notations, we find 

u\N,k,±) = l-2b+^±^V + 4^ , fc = 0,l,2,... (14) 
£ 2 (iV,A;,±) = uj 2 (N,k,±) , k = 1,2,3... (15) 

This is reminiscent to the situation in the elliptic Perelomov model: the elliptic 
model is not solvable for generic values of the parameters determining the period but 
in the limits corresponding to rational or trigonometric potentials the model becomes 
completely solvable. 

The above system is related to the normal mode analysis of the soliton appearing in 
the two dimensional Goldstone-model with periodic boundary condition of the space 
variable [8]. We would like to point out that coupled-channel Hamiltonians which 
posses a multiple algebraisation are rare, very specific and difficult to construct but 
many of them are related to the normal mode analysis of classical solutions (solitons or 
sphaleron) of simple field theories. These algebraic properties can be used as testing 
ground for more ambitious models. 
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